Scaling and dynamics of washboard road 



Anne-Florence Bitbol^'^, Nicolas Taberlet^, Stephen W. Morris^, Jim N. McElwaine^ 

^ Universite de Lyon, Ecole Normale Superieure de Lyon, 
Laboratoire de Physique, 4^ allee dltalie, 69007 Lyon, France 
'^DAMTP, University of Cambridge, Wilberforce Rd., CBS OWA Cambridge, U.K. 
^Department of Physics, University of Toronto, 60 St. George St., Toronto, Ontario, Canada, MBS 1A7 

(Dated: June 19, 2009) 

Granular surfaces subjected to forces due to rolling wheels develop ripples above a critical speed. 
The resulting pattern, known as washboard or corrugated road, is common on dry, unpaved roads. 
We investigated this phenomenon theoretically and experimentally, using laboratory-scale apparatus 
and beds of dry sand. A thick layer of sand on a circular track was forced by a rolling wheel on an 
arm whose weight and moment of inertia could be varied. We compared the ripples made by the 
rolling wheel to those made using a simple inclined plow blade. We investigated the dependence 
of the critical speed on various parameters, and describe a scaling argument which leads to a 
dimensionless ratio, analogous to the hydrodynamic Froude number, which controls the instability. 
This represents the crossover between conservative, dynamic forces and dissipative, static forces. 
Above onset, wheel-driven ripples move in the direction of motion of the wheel, but plow-driven 
ripples move in the reverse direction for a narrow range of Froude numbers. 
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The spontaneous rippling of gravel roads after the pas- 
sage of many vehicles is annoyingly familiar to drivers of 
unpaved roads around the world. Avoiding or mitigat- 
ing the rippling effect, known as washboard or corrugated 
road, is a significant engineering challenge [U El O HI 
m El |7]. We have argued recently |8] that the rippling 
can be regarded as a type of nonlinear pattern-forming 
instability [9] of the flat road [lOl [11]. This nonlinear 
physics point of view is in many ways complementary 
to the engineering one, and leads to new insights. We 
present here an experimental and theoretical study of 
the region near the onset of the instability, which ap- 
pears only above a critical threshold speed. We focus 
on how this critical speed scales with various parame- 
ters. We also survey some of the complex dynamics that 
occurs above onset [10]. 

Washboard road is commonly, and incorrectly, as- 
sumed to require the forced oscillation of the suspen- 
sion of the vehicle [7|. Although a suspension system 
does change the quantitative details, ripple frequencies 
are generally quite far from the free resonant frequencies 
of the suspension [2^ X, , and ripples can appear even with 
no suspension present [8]. Here, we consider a simplified 
experimental system in which a wheel merely returns to 
the road surface under gravity. We also study the even 
simpler case of a non-rotating, angled "plow" blade in 
place of the rolling wheel. In both cases, well above on- 
set, we find that the ripples travel down the road in the 
direction of driving and that the wheel or plow is thrown 
free of the roadbed between ripples. In the case of the 
plow, there also exists a range just above onset where the 
blade remains in contact with the surface and the ripples 
travel in the reverse direction, against the direction of 
driving. 

The rippling instability is quite robust, and has analogs 
in other similar effects on granular surfaces, such as 



wind ITT and water-driven \TT ripple patterns. Other re- 
lated phenomena are, for example, periodic wear patterns 
on railroad tracks [14] , mogul formation on ski slopes [15] , 
skipping stones over water [16] and certain failure modes 
of computer hard disks [TTl [18] . What all these phenom- 
ena have in common is a repeated, nonlinear interaction 
between a moving source of lateral and vertical contact 
stresses acting on a deformable or erodible surface. The 
interaction includes the case when the source of stress, 
which both deforms and responds to the shape of the 
surface, completely leaves the surface so that the stresses 
go to zero. We propose a very general picture of the 
key mechanism of the washboard instability in terms of 
a dimensionless group, analogous to the hydrodynamic 
Froude number, which we show controls the scaling of 
the onset speed. This analysis may be broadly general- 
izable to the diverse situations cited above. 

Previous experimental work and soft sphere molecular 
dynamics simulations [8^ have shown that neither com- 
paction nor size segregation processes in the granular bed 
are essential for the instability, contrary to some theo- 
ries [11,. Engineering studies [1, 2^ 3^ 4^ 5^ 6 have mainly 
concentrated on optimizing the response of a suspension 
system to fully developed ripples and have not typically 
considered the threshold of the instability. Using our ap- 
proach, it may be possible to suppress or delay the onset 
of this threshold by using appropriate control strategies 
based on the nonlinear dynamics of the coupled wheel- 
roadbed system, which behaves in some respects like an 
impact oscillator [T9l [20] . 

This paper is organized as follows; in the section [l| 
we describe the experimental apparatus. This is followed 
in section [ll| by a qualitative discussion of the rippling 
phenomena that are observed for both wheels and plow 
blades. In section |lll| we consider the theory of the basic 
scaling of the onset speed, which we then examine exper- 



2 




FIG. 1: The two versions of the experimental apparatus. The left photo shows the setup in Cambridge, where the table carrying 
the sand roadbed rotates beneath an arm holding a wheel or plow, which is stationary in the lab frame. The right photo shows 
the setup in Lyon, where the roadbed is stationary beneath an arm that rotates. In both experiments the moving parts rotate 
counterclockwise, as seen from above. 



imentally in section |IV| Next, in section [V| we discuss the 
regime of backward traveling ripples and in section VI 
we describe some of the rich variety of ripple states that 
appear above onset. Finally, section VII is a brief con- 
clusion. 



I. EXPERIMENTAL METHOD 

Two versions of the experimental apparatus were used 
in this study, as shown in Fig. [l] The first was essentially 
that of Ref. [8] and consisted of a 1 m diameter rotating 
table, with a maximum rotation speed of 0.8 Hz, holding 
a 200 mm deep layer of rough sand, producing a moving 
roadbed. This passed under an arm which could hold 
either a rolling wheel or a plow blade which was sta- 
tionary in the lab frame. This arrangement has some 
advantages for visualization and also the large moment 
of inertia of the table helps maintain a constant speed. 
In the second version of the apparatus, the arm hold- 
ing the wheel or plow hangs from a rotating axle over a 
1.80 m diameter roadbed which is stationary in the lab 
frame. This version is larger and can reach higher speeds, 
with a maximum axle rotation frequency up to 1 Hz. In 
both versions, we used rough sand with a grain diame- 
ter of 300 ±100 fim. We have previously shown [8^ that 
the shape and size of the grains is unimportant to the 
rippling phenomenon. The density of the grains, which 
was not varied in the present study, is however a crucial 
parameter, as discussed in section [Hi] below. 

Fig. [2] shows the arm with the wheel and plow blade 
configuration in detail. A light, 330 mm long arm pivoted 
around a point of support. An optional counterweight 
could be added to the end opposite the wheel or plow. 
The weight and transverse width of the wheel or plow can 
be varied, and, in the case of the plow, the angle of attack 
a was also adjustable. The wheel used was smaller and 
lighter than that described in Ref. [8] , having a diameter 



of 65 mm and a typical width of 25 mm. The plow blade 
was typically 120 mm wide and inclined at an angle of 
a = 40°. 

The masses of the wheel or plow and counterweight 
could be varied between 0.045 kg and 3 kg. As discussed 
in section [HT] below, these arrangements are dynamically 
equivalent to a physical pendulum, as shown in Fig. 
Since the angle of the arm was small, the physical pen- 
dulum was always nearly horizontal. Thus, the wheel or 
plow behaved as if it were a freely falling object and the 
natural frequency of the pendulum played no role. No 
spring or dashpot was used, so this arrangement consti- 
tutes a minimal suspension having no resonant or natural 
frequencies, and only rather weak damping. 

The wheel rolled freely and no torques were applied 
to it other than that produced by contact with the sur- 
face. The plow was simply dragged over the surface. The 
speeds were inferred from measurements of the rotation 
rate of the table and spanned the range 0-6 ms~^. A 
magnetic angle sensor [2l was attached to the pivot sup- 
porting the arm, and its angle 6 was digitized and used 
to find the vertical position z of the wheel or plow blade. 
z was measured with a vertical accuracy of 0.2 mm at 
a rate of 2000 Hz. The position data was acquired con- 
tinuously along with data from an optical encoder that 
generated a pulse each table rotation. 

The dynamics of the ripple pattern evolved rather 
slowly, especially near the critical speed Vc for the on- 
set of the instability. It was thus necessary to increase 
and decrease the speed over periods of many hours to 
achieve quasi-static, steady state patterns. In some runs, 
the speed was ramped up quickly from a standing start, 
beginning with a flat sand bed. This protocol is analo- 
gous to making a rapid "quench" into the rippled state. 
In other runs, a small perturbation was made in the bed 
to act as a nucleation site for ripples. 
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FIG. 2: Schematic views of (a) the wheel and (b) plow blade, 
(c) shows the geometry of the equivalent physical pendulum. 
The plow is rigidly fixed to the arm and makes an angle of 
attack a with respect to the horizontal. 



II. QUALITATIVE RESULTS 

In this section, we discuss some of the phenomena ob- 
served, before focusing on the scaling of the critical speed 
for the onset of rippling in section [TV} Initially, we tried 
preparing the roadbed by running a comb through it to 
make a loose, fiat, granular assembly. However, if the 
experiment is run continuously over a period of around a 
week 100 000 rotations), the bed gradually compacts 
and the thresholds and amplitudes slowly change. We 
therefore used a different protocol where the sand was 
disturbed as little a possible. To prepare a flat bed, the 
experiment was run at a speed below the onset for many 
rotations. The results were then repeatable. 

Ripples driven by either the wheel or the plow only 
appear above a well-defined critical speed Below this 
speed, perturbations made in the roadbed are observed to 
be smoothed away by the action of the wheel or plow. For 
speeds near the growth rate of perturbations is nearly 
zero, and it can take many thousands of rotations of the 
experiment before ripples build up. The quasi-steady 
state ripples which develop after the initial transient have 
a very small amplitude just above For all cases, speeds 
well above Vc result in regimes where the wheel or plow 
is thrown completely off the roadbed near the crest of a 
ripple, before landing again on the front face of the next 
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FIG. 3: (Color online) Typical space-time plots showing the 
development of ripples from a flat bed, with red indicating a 
higher arm position and blue a lower. Both parts show 1000 
rotations and have the same scale. The driving direction is to 
the right, a) A free wheel just above the critical velocity, b) 
The plow, with v below the critical velocity for the first 300 
rotations and well above critical for the next 700 rotations. 



ripple. The resulting large amplitude ripples are strongly 
asymmetric, with a gradual rise on their front faces, and 
an abrupt slip-face on the back side. Nearer to the 
ripples are more symmetrical. 

The plow is observed to remain in contact with the sur- 
face at all times near Vc- The ripples move in the direction 
of driving for wheels, but may move in the reverse direc- 
tion for plows, if V is close to critical. At larger speeds, 
ripples move in the direction of driving in both cases. We 
discuss this curious reversing behavior in section [V| be- 
low. Fig. [3] shows some representative space-time plots of 
ripple evolution after a quench into the strongly unsta- 
ble regime. During the transient, numerous ripple cre- 
ation and annihilation events are observed, resulting in 
adjustments of the pattern which tend to propagate in 
the direction of driving. Fig. [4^ shows a quench close 
to the threshold with the plow showing backward rip- 
ples. Eventually, a quasi-steady state of traveling ripples 
is established. For sufficiently deep quenches, the pattern 
can display multi-modal behavior or persistent time de- 
pendence, as discussed in section [Vl The final, saturated 
amplitude of the ripples increases with the driving speed 
near the threshold, but shows non-monotonic behaviour 
far above. 

It is easy to observe experimentally that Vc depends 
on the mass of the wheel or plow, and the moment of 
inertia of the arm. By changing the granular material, 
it was apparent [8 that only the density of the roadbed 
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FIG. 4: (Color online) Parts a and b show different segments 
of a ten day experiment during which the speed was decreased 
in steps of 0.0021 ms~^ every 500 rotations. The driving di- 
rection is to the right. In part a, the speed is increasing 
from 0.684 ms"^ to 0.691 ms"^ (Fr 2.04-2.09), triggering a 
transition from flat bed to backward traveling ripples. In 
b, The speed is decreasing from 0.694 ms~^ to 0.687 ms~^ 
(Fr 2.10-2.06), triggering a transition from 21 to 20 ripples. 
Here Mq 0.0785 kg, w 0.120 m, p 1201 kg m^ and the 
angle of attack a was 36.8 ° 



material matters, and not the size or shape of the grains. 
These dependencies can be accounted for by the scaling 
analysis discussed in the next section. 



III. SCALING THEORY 

In this section, we consider the dynamics of the arm 
assembly and its interaction with the sand bed. We seek 
to identify the essential physical parameters that gov- 
ern the scaling of the critical onset speed Vc- To sim- 
plify the problem, we consider a strictly two-dimensional 
model in which the transverse width of the wheel or 
plow is assumed to be large compared to its penetra- 
tion into the bed. This assumption was exactly true for 
the soft sphere molecular dynamics simulations discussed 
in Ref. [8], which qualitatively reproduced all of the es- 
sential features of the experiment. The validity of this 
simplifying assumption for the real experiment will be 
discussed in section [TVl 

The geometry of the arm and plow assembly is shown 
in Fig. The position of the assembly is completely 
specified by three points: the contact point with the bed, 
the pivot position and the location of the center of mass. 



These three positions define two intrinsic lengths, the dis- 
tance from the pivot point to the center of mass ^, and 
the distance from the pivot point to the contact point 
with the bed L. The angle (j) between these lines is con- 
stant. The angle a — between the plow and the line to 
the pivot point is also fixed, whereas 6 will vary as the 
bed deforms. 

Dynamically, the system can be described as a forced 
physical pendulum of length ^, mass m and moment of 
inertia about the center of mass /. The latter two quan- 
tities define another length ^Jl|m. There is one addi- 
tional length, the transverse width of the plow or wheel 
w. Thus, to describe the whole assembly requires three 
angles, a, Q and (which are of course non-dimensional), 
four lengths ^, L, ^Jl jm and the mass m and the ac- 
celeration due to gravity g. The properties of the bed 
involve only one additional parameter, the sand density 
p, as discussed below. 

The gravitational force mg acts through the center of 
mass and exerts a torque mgicos{0 — (j)) about the pivot 
point. The force of the bed on the arm is assumed to act 
entirely at the contact point. It exerts a torque about 
the pivot point of —L{N cos + Tsin^), where N is the 
vertical (normal) force from the bed and T the horizontal 
(tangential) force. Since the moment of inertia about the 
pivot point is /+m^^, conservation of angular momentum 
gives 

(/ + mf)0 = mg£ cos(0 - 6>) - L(7V cos (9 + T sin 0). (1) 

The vertical height of the contact point relative to the 
pivot point is z = —LsinO. To determine how the onset 
speed Vc scales, we consider small displacements about 
a flat bed, so that is small. This is sufficient, since 
the transition to rippling occurs at zero amplitude. We 
treat as constant except in the term where z = —LO. 
After dividing Eqn.[l]by L and substituting for z we get 
approximately 



Mi 
cos 6 



z = {NcosO ^TsinO) - Mg gcos{0-(j)), (2) 



where Mj = (/ + mP)/LF' is the effective inertial mass 
and Mq = ml/L is the effective gravitational mass. 

Mj and Mq can be changed independently by adjust- 
ing ^, L and m. They are equal if all the mass is located 
at the end of the arm. If the angular dependencies are ab- 
sorbed into the definitions, that is with M'j = Mi /cos 6^ 
M^ = cos{0 - (}))Mg, and N' = N cosO ^ TsmO, then 
for small displacements {i.e. small changes in 0) we have 
the equation of motion 



M'jz = 



(3) 



which is simply that of a falling mass interacting with a 
surface via a normal force N' . This equation holds for 
small 6>, but for larger ^, the coupling of the resistive force 
T to the vertical motion becomes important. Our exper- 
iments concentrated on small ^, but in subsection IV C[ 
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we discuss the results of experiments where this condition 
is relaxed. 

The mechanics of the granular road bed are not so 
straightforward to model as those of the arm. We aim for 
a minimal description that leads to the correct scaling of 
the fully coupled problem. Much more detailed analyzes 
of the rheology of the granular material would be required 
for a very realistic description [3j. 

We consider the case where the depth of penetration 
into the bed is small compared to the width of the object, 
so that we may assume that the induced flow in the sand 
will be essentially two-dimensional, and that edge effects 
are small. This assumption is probably violated in our 
experiments with the narrowest wheels and plows. 

The sand grains were dry and sufficiently large so that 
they were non-cohesive. Previous work [8 has shown that 
changing the grain size does not change the wavelength, 
amplitude or critical velocity of the ripples. Therefore, 
the granular medium can be described by only its three- 
dimensional density p and possibly by some other non- 
dimensional parameters. The shallowness assumption 
and two-dimensionality means that the density p of the 
sand and the width w will only occur in the form of a 
two-dimensional density pw^ and that the width should 
have no further role. It may be possible to allow for some 
three dimensional effects by allowing a virtual origin cor- 
rection in the width, though we have not tested this. 

To arrive at a plausible scaling, we examine the balance 
of forces for a plow or wheel moving over and penetrating 
the sand surface. Forces between the bed and the plow 
or wheel arise in two distinct ways which scale differently 
with V. One of these is almost entirely dissipative and 
the other one is almost entirely conservative. We postu- 
late that the onset of washboard ripples occurs when the 
relative sizes of these forces cross over as a function of v. 
We consider the geometrically simpler case of the plow 
first. 

Static, buoyancy or frictional forces must be propor- 
tional to Apwg^ where A is the displaced area of granular 
material. For a plow penetrating a distance h = zq — 
where zq is the free surface position, they will be pro- 
portional to h^pwg. These forces will be almost entirely 
dissipative. On the other hand, dynamic forces, those 
arising from the inertia of the sand, will be proportional 
to hpwv^. These forces are mostly conservative. We ne- 
glect the vertical velocity h as small compared to v. Ig- 
noring constants and angular dependence, we separately 
balance these two penetration forces against the effec- 
tive weight of the arm Meg- This gives the static and 
dynamic equilibrium penetration depths 
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hd 




v'^pw 



(4) 
(5) 



As V increases, the dynamic forces increase and the dy- 
namic penetration depth hd needed to support the weight 



decreases, whereas the static forces are unaffected. Since 
the static forces are almost completely dissipative and 
the dynamic forces almost completely conservative, the 
switch between the two that occurs as v increases ex- 
plains the transition from a highly dissipative furrow- 
plowing mode to one in which the plow spends some 
of the time in non-dissipative, ballistic flight. For low 
speeds the normal forces between the plow and the 
sand do only dissipated work against friction, while for 
larger the normal force is sufficient to throw the plow 
off the surface and into free flight. Ripples are then 
formed and maintained by the periodic dissipative col- 
lisions with the surface. The onset of washboard ripples 
in this picture is thus somewhat analogous to the onset of 
skipping for a stone projected over a water surface [16 . A 
key difference, of course, is that the granular ripple pat- 
tern is built up by many passes of the plow, each of which 
contributes only a small, but persistent, deformation of 
the surface. The important non-linearity in the problem 
is the highly nonlinear nature of the normal force, which 
abruptly goes to zero when the plow leaves the surface, 
as is the case for impact oscillators [191 EH] • 

The dimensionless ratio of the two penetration depths. 
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is analogous to the Froude number used in similar hy- 
drodynamic arguments [16 . We propose [8 that this 
dimensionless group controls the onset of the instability, 
with ripples setting in above a critical, threshold value of 
Fr = Frc. Our scaling argument does not furnish any es- 
timate of the value of Frc, beyond the general expectation 
that Frc ^ ^(1)- We establish experimental values of Frc 
in section [TVl below. We will also show cases where the 
plow does not actually leave the surface, but still forms 
backward traveling ripples for Fr > Frc • 

We have physically motivated the identification of the 
Froude number in terms of a crossover between dissi- 
pative and conservative forces [22 , but it can also be 
deduced from straightforward dimensional analysis. If 
the only relevant dimensional quantities of the arm that 
determine the transition are its effective gravitational 
mass Mq and its width and if we assume the two- 
dimensional model such that w and p can only appear as 
the product pw^ then Fr is the only dimensionless group 
that can be formed from a combination of M^, pw^ v and 
9- 

The critical Froude number Frc will also depend on 0, 
a, (j) and possibly on other dimensionless properties of 
the bed. If these are held constant, Eqn. [6] implies that 
the critical speed Vc is given by 



\pw J 



1/4 



(7) 



In subsection |IV A[ we will show experimentally that Vc, 



for the plow, does indeed exhibit an Mq^ scaling. In 



6 



Sections |IVB| and |IVC[ we consider the dependence of 
Vc on a and 6. 

The wheel is more comphcated than the plow in some 
ways as there is an additional length scale, the wheel 
radius R. The displaced volume for a cylinder length w 
which penetrates to depth h is 



w 



cos-\l - h/R) -{R- h)^2Rh-h^ 



= -whVWh 



0{h/R) 



(8) 



Thus the scaling is now A ex h^^'^ instead of A ex /i^, so 
that the static depth is modified. For the wheel 
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(9) 



The dynamic depth is the same thus the scaling of the 
Froude number is changed to 
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(10) 



Thus, we predict the following relation for the critical 
speed 



..wheel 
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(11) 



where numerical coefficients have been ignored. 

To experimentally determine the relevant parameters 
of the model, it suffices to determine the three constant 
quantities that appear in Eqn.[2j namely M/, Mq and (j). 
We supported the arm on a balance, so that the bottom 
of the plow was exactly level with the pivot point and 
^ = 0. Then the force measured by the balance was N = 
Mq gcosc/). The angle (j) was determined by suspending 
the arm from the pivot point and measuring 0. Then 
Mg gcos{e - ^) = and = ±7r/2 + 6. Eqn. [2]for the 
suspended arm is simply the equation of a free physical 
pendulum. 



MiL = Mg gcos{(j)-0), 



(12) 



which, for smal l amplitude s, oscillates with angular fre- 
quency uj — ^Mg g/MiL. The length L was directly 
measured with a ruler, while co was measured by trig- 
gering small oscillations, and then measuring using the 
magnetic sensor. In this way, all the parameters in Eqn.|2] 
were determined. All these quantities vary slightly with 
the angle of attack a of the plow, and so were measured 
for a range of a. The plows of different width were care- 
fully chosen so as not to alter any of the properties. 



IV. SCALING OF THE THRESHOLD SPEED 

To test the scaling theory outlined in the previous 
section, we measured the critical speed Vc for a vari- 
ety of plow and wheel configurations. The transition 



speed is most accurately determined by running the sys- 
tem at a constant speed above Vc until a steady state 
has developed. The speed is then decreased by steps of 
0.0016 ms~^ and kept at each speed for 1000 rotations. 
This is more accurate and quicker than starting with an 
initial perturbation or fiat bed and observing the sub- 
sequent behaviour, as a very large number of rotations 
are necessary for the ripples to develop near Vc. We also 
recorded the long-term dynamical states for a range of 
speeds from below to well above onset. 

As discussed in the following subsections, we added 
various masses and counterweights to the arm, which has 
the effect of varying Mj and M^, as well as varying the 
width of the plow its angle of attack a and the mean 
support angle 0. 



A. Variation with gravitational mass Mg 
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FIG. 5: (Color online) The domain of existence of the rip- 
ple pattern for the plow, as a function of the speed v and 
the gravitational mass Mg- The solid line corresponds to a 
constant Froude number Fr = 2.58. The inset shows v'^ / g vs. 
[Mg / pwY^'^ ^ a scaling for which Eqn.[6]predicts the boundary 
to be a straight line. 

Fig. [5] shows the variation of Vc with the gravitational 
mass Mg for the plow. The frontier between the ffat bed 
and the appearance of washboard ripples is well described 
by an Mq^ scaling, corresponding to a critical Froude 
number of Fr = 2.58. This result clearly validates the 
theory presented in the previous section, and indicates 
that the instability is controlled by Fr. It is interesting to 
note that the first patterns encountered above onset are 
backward traveling ripples for which the plow does not 
lift off the surface. These are discussed in more detail in 
section IVl below. 

The scaling theory is less successful for the case of the 
rolling wheel; Fig. [6] shows Vc vs. Mg for the wheel. To 
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FIG. 7: The dependence of the critical velocity Vc on the 
width w of the plow. The Froude numbers of the three points 
shown were 1.84, 1.86 and 1.99, and hence the critical Froude 
number Frc remained approximately constant while w was 
varied. The solid line shows a satisfactory fit to the w'^^^ 
scaling expected from Eqn. jT] 



of the granular surface, presumably scale with additional 
powers of R and perhaps also depend on the moment of 
inertia of the wheel. In any case, Fig. |6] shows that the 
gravitational mass dependence of Vc is extremely weak 
for wheels, rolling or otherwise. 



FIG. 6: (Color online) The critical speed Vc for the onset of 
the ripple pattern for wheels. Part (a) shows Vc vs. Mq for 
"blocked", or non-rolling wheels, while (b) shows the same 
data for freely rolling wheels. The solid line corresponds to 
a constant Froude number, according to Eqn. |11[ i.e. to an 
Mq^ power law dependence. It is evident that the predicted 
scaling is not obtained for rolling wheels. 



isolate the effects of rolling from the purely geometrical 
dependence, we measured Vc both for wheels that were 
free to roll and for otherwise identical wheels that were 
"blocked" — i.e. prevented from rotating. A blocked 
wheel behaves much like a differently shaped plow. 

Fig. [g^ shows that the predicted Mq^ scaling is con- 
sistent with the blocked wheel data. However, exper- 
imental limitations only allowed us to vary Vc over a 
range of 20%, so the correctness of this scaling is not 
conclusively demonstrated. It is clear however that this 
scaling fails for the case of a rolling wheel, as shown in 
Fig. [gJd. This is perhaps not surprising since arguments 
presented in section [TTT| were based only on the geometry 
of the wheel and completely ignored its rotation. Evi- 
dently, the rolling component of the motion affects the 
surface penetration properties of the wheel, invalidating 
the scaling based only on geometry. The contact forces 
for a rolling wheel, and their effect on the deformation 



B. Variation with plow width w 



The scaling argument presented in section [IV| de- 
pended on the assumption that the motion was essen- 
tially two-dimensional so that only the two-dimensional 
density pw mattered. This is probably the most drastic 
approximation in the scaling theory. Figure [7] shows how 
Vc scales with w for the plow, for constant a and Mq- 
The data is well fit by a w~'^^^ dependence, as predicted 
by Eqn. [Tj and we do find that the critical value of Fr 
remains constant even as w is varied by a factor of four. 
The w scaling can only hold for a relatively small range 
of plow widths, however. For very narrow plows, w ^ h 
and the two-dimensional scaling will break down. On the 
other hand, for very wide plows w becomes a significant 
fraction of the radius of the circular bed and thus there 
is a significant variation of v across the width. Figure [7| 
spans the intermediate range of w where a reasonable 
agreement is obtained. 



C. Variation with the angles a and 

The onset of the washboard ripples depended strongly 
on the angle of attack of the plow a, and on the arm 
angle 6>, as defined in Fig. [2j a could be varied, while 
keeping all other parameters constant. depended on the 
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FIG. 8: The dependence of the critical Froude number on 
the angle of attack of the plow. The angle a is measured 
with respect to the horizontal (see Fig. [2]). At low angles, the 
transition was very sharp and the error bars lie within the 
symbols. 



FIG. 9: The dependence of the critical Froude number on the 
mean angle that varies with the height of the support. 



forces, we would guess that 



height of the supporting structure, and could be adjusted 
independently of a. Due to experimental limitations, 9 
could only be varied over a range of about 15^. The 
two angles were in fact somewhat interdependent, as the 
mean height of the bed could slowly change over time 
due to compaction or to sand moving to other parts of 
the table. In such cases, a and both drifted upward by 
a few degrees over the course of an experiment. 

The simple scaling theory of Eqn. [7| does not explic- 
itly predict the a and dependence of the critical speed 
which would be expected to show up as a variation 
of the critical Froude number Frc. Fig. |8] shows how 
Frc depended on a. We found that and hence Frc, 
strongly decreased as a increased. This can be explained 
by the observation that the transport of sand along the 
bed increases when a is increased. It is not a surprise, 
therefore, that the velocity needed to trigger the wash- 
board instability should decrease with increasing a. The 
case of a ~ 90°, a nearly vertical plow, is somewhat pecu- 
liar. At this point, the transition becomes poorly defined 
and the amplitude noise near onset greatly increases. At 
a ^ 90°, it appears that other parameters such as the 
thickness of the plow may play an important role. 

The systematic interdependence of a and makes it 
quite difficult to exactly repeat the experimental condi- 
tions. The error bars in Fig. [S] are the precision that the 
transition can be identified within a particular experi- 
mental run rather than the somewhat larger error bars 
that would be obtained from trying to achieve exactly 
the same 6 over many runs. 

The critical speed Vc also shows a strong dependence 
on 0. Fig. |9] shows how the critical Froude number scales 
with variations in 6>, with a = 45 ± 3°. If we directly 
apply Eqn. [T] and make a naive model of the friction 
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where /i is the unknown ratio between the normal force 
N and the tangential force T. For 0(1) values of /i and 
small ^, this predicts that the critical Froude number 
should decrease weakly with increasing 0. This is the 
opposite of what is observed. Clearly, the dependence 
of the critical Froude number is not so easily explained 
and a more detailed study is required to account for it. 



V. BACKWARD-TRAVELING RIPPLES 

In all experiments using wheels, including our previ- 
ous studies [8 , the ripples that appear just above onset 
were observed to travel in the direction of driving. This 
agrees with most, but not all, observations in previous 
engineering studies [1 . However, ripples produced by a 
plow just above onset always travel in the reverse direc- 
tion, i.e. opposite to the direction of driving, as shown 
in Fig. |4] The backward travel reverses at higher driving 
speeds, but there is no smooth transition from backward 
to forward ripples; there is no speed where the ripples 
are stationary. If the speed is rapidly stepped from be- 
low to well above as in Fig. [3]3, the threshold speed 
above which backward ripples are observed in the steady 
state appears to scale in the same way as Vc, as shown 
Backward traveling ripples only exist for a 



in Fig. 10 



narrow range of parameter values just above onset. 

Another distinguishing feature of backward traveling 
ripples is that ripple merger events, which are common 
near onset for forward traveling ripples, generally do not 
occur. This is because the backward ripples appear at the 
same wavelength as in their amplitude-saturated state. 
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whereas there is a strong relationship between the wave- 
length and amplitude for the forward ripples. This is 
related to another difference: that the plow does not 
leave the granular surface during backward ripple motion, 
while such ballistic motion is always present for forward 
traveling ripples using the plow. 




^0 100 200 300 400 500 600 



on the upward face of the ripple than on the downward 
one. The upward face is therefore less eroded, while the 
sand that is preferentially eroded from the downward face 
is deposited on the upward face of the following ripple. 
Thus, the net motion of the ripple is backward, while the 
grains are still transported forward. A similar mecha- 
nism has been proposed for the uphill motion of moguls 
on ski slopes [23 . 

Conversely, in the case of forward traveling ripples, 
inertial effects dominate and the plow leaves the surface 
and goes into ballistic motion near the crest of the ripple. 
The erosion is therefore all on the upward face, while the 
downward one forms a slip face which is added to by 
small avalanches consisting of material transported from 
the upward face. Thus, the ripples in the ballistic regime 
move forward, in the driving direction. 



FIG. 10: (Color online) The domain of existence and tran- 
sitions for forward and backward traveling ripples using the 
plow. The data and upper solid line shows the transition be- 
tween backward and forward traveling ripples at Fr = 3.80. 
The lower solid line shows the onset of backward ripples from 
the flat bed at Fr = 2.58, as in Fig. [5] 




FIG. 11: Schematic mechanism for the formation of backward 
traveling ripples. The arrow shows the direction of motion of 
the plow. If the rate of erosion by the plow increases with 
the angle (3 with the surface, then material is preferentially 
removed from the downward face of the ripple and deposited 
on the upward face of the following one, giving rise to a net 
backward ripple motion. During this process, the plow does 
not leave the surface. 

The formation of backward traveling ripples is funda- 
mentally different from the inertial mechanism of for- 
ward traveling ripples. Ripples can move backward even 
though the sand is only transported forward. The sur- 
face flux of grains in the driving direction depends not 
only on the normal force, but also on the geometry of the 
contact region. It is plausible to suppose that the sur- 
face flux increases with the angle /3 between the plow and 
the granular surface. Initially, /3 = a, the fixed angle of 
attack of the plow, but as the ripples develop, P is mod- 
ulated by the ripples as the plow remains in contact with 



VI. MORE COMPLEX STATES 



In the preceding sections, we have focused on the scal- 
ing of the critical speed Vc for the onset of the ripple 
pattern. While a complete study of the states above on- 
set is beyond the scope of this paper, it is interesting 
to survey some of the complex states that are encoun- 
tered above Vc. Fig. [12] shows a selection of such fully 
developed ripple states. 

Other than the flat bed, the simplest states are pure 
traveling ripples, which are well described by single 
Fourier modes. For wheels, these always travel forward, 
while for the plow they travel backward near onset, as in 



Fig. [T2p, and forward for larger speeds, as in Fig. [T2p. 
As the driving speed v of the plow is increased, the back- 
ward traveling ripple speed decreases, but is never zero. 
Increasing v results in a low-dimensional chaotic state, as 
in Fig. [l2]3, in which several Fourier modes are excited. 
Further increasing we enter a regime of steady forward 
traveling ripples, but here the plow spends much of the 
time in a free, ballistic trajectory between impacts. If 
these impacts are not periodic around the table, there is 
some jitter between the phase of the ripples and that of 
the table; this is visible in Fig. 12z. At still higher speeds 



of the plow, a second chaotic regime is reached, as shown 
in Fig. [I2li. 



The phenomenology for wheels differs in detail, but is 
similarly complex. For example. Fig. 12 3 and f show 



the surface, as shown in Fig. 11 The result of this mod- 
ulation is that and hence the surface flux, is smaller 



wheel-driven ripple states consisting of two significant 
Fourier modes which are not simply related. The two 
ripple modes propagate with different speeds, producing 
a beating effect which persists indefinitely. All of these 
complex states for plows and wheels are reminiscent of 
the multiple dynamical states of impact oscillators |20] 
in their strongly nonlinear regime. 
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FIG. 12: The zoo of states for the plow and wheel. Each figure contains 1000 revolutions at constant speed. The range of ripple 
heights is normalized for each figure with red representing higher a) Plow at 0.737 ms~^. b) Plow at 0.827 ms~^. c) Plow at 
0.964 ms"^ d) Plow at 1.001 ms"^ e) Wheel at 1.166 ms"\ f) Wheel at 1.325 ms"^ 



VII. CONCLUSION 

We have studied the rippling instability of a flat gran- 
ular surface under the action of a moving source of stress. 
We considered both the classic case of a rolling wheel pro- 
ducing "washboard road" and the simpler case of a flat 
plow blade. In both cases, we focused on the region near 
the onset speed Vc. We developed a scaling theory for Vc 
based on the idea that the onset of the instability rep- 
resents a crossover between conservative, dynamic forces 
proportional to v'^ and dissipative, static forces which 
are independent of v. We found that Vc did scale with 
the gravitational mass as Mq as predicted for the plow 
and for non-rolling wheels. These have different scalings 
due to differences in their shape. However, the predicted 
mass scaling was not found for rolling wheels, suggesting 
that a more detailed theory of the surface penetration 
of a rolling wheel is required. We then considered the 
dependence of Vc on various other geometric parameters. 

Finally, we surveyed some of the complex ripple dy- 
namics that are observed above onset, which include both 
forward and backward traveling ripples, chaotic states 
and states with traveling defects. Clearly, understanding 



the above-onset regime of finite amplitude ripples and 
their interactions will require a fully nonlinear theory of 
the instability. Such a theory would need to go well be- 
yond the simple scaling of the onset speed Vc that we 
have discussed here. 



Throughout this work, we have used the angular posi- 
tion of the arm to trace the motion of the wheel or plow. 
To completely close the feedback loop between the mov- 
ing source of stress and the state of the granular bed, 
we would also need to measure the shape of the granular 
surface and the rotational state of the wheel. This would 
allow us to determine the full acceleration of the wheel or 
plow, and hence the forces on it, as well as the dynamic 
response of the granular bed. The bed shape could be 
determined, for example, by a scanning laser profilome- 
ter [24 . Another, somewhat complementary avenue to 
the same kind of information is soft sphere molecular dy- 
namics simulation [8] in 2D or 3D. With this data, a 
realistic dynamical model of the wheel-bed interaction 
could be developed, leading to a more complete contin- 
uum theory of the washboard instability. 
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